Heterogeneous structure of granular aggregates with capillary interactions 
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We investigate the spatial structure of cohesive granular matter with spheres floating at an air- 
liquid interface that form disordered close packings with pores in between. The interface is slowly 
lowered in a conical container to uniformly compress and study the system as a function of area 
fraction </>. We find that the free area distributions associated with Voronoi cells show significant 
exponential tails indicating greater heterogeneity compared with random distributions at low (f) with 
a crossover towards a F-distribution as <j) is increased. Further, we find significant short range order 
as measured by the radial correlation function and the orientational order parameter even at low 
and intermediate (j>, which is absent when particles interact only sterically. 

PACS numbers: 45.70.-n, 68.03.Cd 



The effect of cohesive interactions on the structure of 
granular matter is an important question which has bear- 
ing on its strength and stability. Further, aggregation of 
cohesive particles is a broad physical phenomenon occur- 
ring among others in colloidal systems, dust, and soot 
due to electrostatic interactions, inter-galactic dust un- 
der gravitation, and in powders due to humidity. Recent 
numerical studies have shown that the connectivity and 
rigidity transition for cohesive frictionless particles oc- 
cur at distinct and lower concentration compared with 
non-cohesive particles [1] , and that cohesion and friction 
leads to more porous structures [21 13 ■ While there have 
been some experiments which have examined the effect 
of cohesion on packing of wet granular matter at the par- 
ticle level in the dense regime [3] , there have been no ex- 
periments which have examined structure over a broad 
range of packing fraction. Complex structures have been 
observed with magnetic beads, but the interactions there 
are dipolar [S]. 

Here, we use capillary attraction between spheres float- 
ing at the air-liquid interface [SI [7] as a model granular 
system with isotropic cohesive interactions. Because of 
the exponential nature of attraction between floating par- 
ticles with distance, the interaction is essentially short 
range for millimeter sized particles. Further, a minimum 
force is required to overcome the rolling or sliding fric- 
tion because of the cohesive normal force acting on par- 
ticles in contact. Thermal excitation is insignificant at 
these scales, and therefore floating spheres can be consid- 
ered as a model system to understand the effect of par- 
ticle level interactions on the heterogeneity of cohesive 
granular aggregates without the complexity introduced 
by gravitational gradients. Recently, such a system was 
used to study the granular character of particle rafts by 
compressing them at the boundaries using a Langmuir 
trough [8 , but neither the area fraction nor the particle 
level heterogeneity was examined. 

A novelty of our apparatus is that we can uniformly 
change the particle concentration at the air-liquid inter- 
face by varying the height of a highly viscous liquid inside 




FIG. 1. (a) Schematic diagram of the experimental appara- 
tus. The liquid level inside the funnel is decreased to increase 
the area fraction cj> of the floating spheres, (b) The meniscus 
around floating spheres leads to attraction, (c) Mean num- 
ber of individual aggregates versus (j). Only one aggregate 
remains on average above 4>c ~ 0.44. Sample images of float- 
ing spheres {N = 688) observed at <^ = 0.33 (d), = 0.50 (e), 
and (j) = 0.65 (f), and corresponding Voronoi cells inside the 
dashed circle (g), (h), and (i). Cells corresponding to tightly 
packed regions using criteria d„/d < 0.1 and 0.1 < d^/d < 0.5 
are identified with dark grey/red and grey/green markers, 
where dv is the distance of the particle from the center of 
the polygonal Voronoi cell. 



a conical container. Thus, any two material points on the 
interface approach each other proportional to distance of 
separation as the liquid level is lowered, rather like two 
points on the surface of an inflated balloon from which 
air is released. Using this trick, we compress and vary 
density of the entire system homogeneously. By imaging 
and finding the positions of the particles, we calculate 
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statistical measures including Voronoi cell area distribu- 
tion, and the radial and angular correlations to charac- 
terize the heterogeneity and understand the interplay of 
cohesion and steric interactions. 

A schematic of the experimental apparatus is shown 
in Fig. [ija) and consists of a conical glass funnel filled 
with a 90 — 10 % glycerol-water liquid mixture by mass. 
The hquid has a density pi = 1.235 x lO'^kg-m"'^, dy- 
namic viscosity rj = 0.26 Pa-s, and surface tension 7 — 
6.5 X 10~^N-m^^. The particles used in our experi- 
ments consist of white polyethylene spheres of diame- 
ter d = 3.175mm and density pp = 0.95 x 10^kg-m~^. 
The funnel is treated chemically to create a hydrophobic 
boundary conditions which repels the floating particles to 
about 2 cm from the boundaries. Images of the floating 
spheres are recorded with a 1392 x 1040 pixel CCD cam- 
era, and centers located to within 0.1 d using standard 
particle tracking software. 

The experiments are initialized by first placing N 
spheres at random locations, which are then observed 
to rearrange over a time scale of a few seconds and 
form small aggregates that remain in place for hours [9]. 
This can be explained by the fact that the force be- 
tween spheres separated by a distance I (see Fig. [T|b) 
for small deformations is given by [7] F{1) cx 
— {l/Lc)~^ exp(— Z/Lc) for l/L^ » 1, where the capil- 
lary length Lc = \/^/ Pig ~ 2.3mm < d. By balancing 
capillary and drag forces one finds that time over which 
spheres approach each other increases roughly exponen- 
tially with distance. Next, after a waiting time of 600 s, 
the liquid level is slowly lowered at a constant rate using 
a syringe as indicated in Fig. [TJa) over a time of 300 s [TU] . 
Over 600 s the range of attraction is about 4d, and during 
the compression, the motion for isolated spheres due to 
capillary forces occurs only over roughly Q.2d. Because 
of the symmetry of the conical funnel, the radial velocity 
of a liquid element at the air-liquid interface away from 
the boundaries is approximately given by v{r) ~ —Co 1^, 
where Co is measured to be 3 x 10"'^ s^^ using a tracer 
technique and is the radial vector to the center axis of 
the funnel. Such a fiow results in isotropic homogeneous 
compression of particles floating on the surface, because 
the strain is independent of position. Further, inertial ef- 
fects are negligible because the maximum Reynolds num- 
ber calculated using the particle size is approximately 
5 X 10-3. 

Examples of aggregates observed are shown in Fig.[l]^d- 
e-f) at various stages of compression. The corresponding 
area fraction (j) is obtained as the fraction of area occu- 
pied by the spheres inside a circle (yellow dashed line) 
whose center is the mean position of spheres and radius 
is chosen as ^3/2 times the radius of gyration in or- 
der to avoid any boundary effects. Voronoi tessellation is 
computed using Qhull programs in MATLAB in the bulk 
indicated by the dashed circle indicated in Fig. [l]jd-e-f). 
No particular variation in size with distance from center 
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FIG. 2. Probability distribution function (PDF) of the free 
area A normalized by the median value of A is observed to 
deviate from a T distribution [v — 3.6) postulated for random 
distribution of points. The deviation is strongest at low <j> 
with significant probability of finding large cells, (b-c-d) The 
tail of the free area PDF with all the cells can be fitted by 
an exponential function with a decay parameter a = 0.484 
(b), a — 1.11 (c), and a = 2.82 (d). When Voronoi cells are 
selected with dv/d < 0.5 the exponential tail decreases, and 
the PDF approaches a F distribution — 3.6) for dv/d < 0.1. 



can be noted supporting our claim that the compression 
is uniform [10]. 

The aggregates initially appear to phase separate into 
regions with high density of particles with voids or pores 
in between. As the system is compressed, the number of 
aggregates decreases as shown in Fig. [T](c) till a thresh- 
old is reached where only one connected aggregate re- 
mains. This connectivity threshold is found to occur at 
(j)c = 0.44 ± 0.02 by averaging over several experiments. 
Upon further compression, void spaces decrease in size, 
until a critical area fraction (pj — 0.715 ±0.006 is reached 
where the aggregates appear to buckle, and the system is 
no longer two-dimensional. This buckling transition has 
been studied previously [S], and may be considered as 
the point where jamming occurs for this system because 
the system can no longer accommodate in plane deforma- 
tion upon further compression. With this assumption our 
data provide evidence that the connectivity transition 
occurs at a distinct and lower 0^ compared with [T], 
while the values themselves are not universal and depend 
on the strength of interaction. This is in contrast with a 
non-cohesive system where these points coincide [T]. Fi- 
nally, it is noteworthy that is significantly lower than 
observed with frictional and cohesionless disks [11 , where 
(t)j - 0.84. 

In order to examine the heterogeneity of the packing, 
we analyze the size distribution of the Voronoi cell . 
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In case of random uncorrelated points, the distribution 
of the area of the Voronoi ceU has been argued to be 
a r distribution [13]: 



P(A„) ^ A^-' exp 



< A, > 



where, (Ay) is the mean area and the shape parame- 
ter v ^ 3.6 from geometric arguments. Now for hard 
disks, a free area A is introduced, which is given by 



A — Ay 



^,2 



d , the minimum possible ceU area cor- 



responding to a hexagonal packing. Numerical simula- 
tions [11] with frictionless hard disks have shown that 
the free area distributions are consistent with a F dis- 
tribution, and experiments with steel cylinders |15j have 
shown V between 3 and 4 depending on the vibration and 
the roughness of cylinder. From Fig. [TJd-f-h) large varia- 
tion in the sizes of the cells are apparent at low (f>, but less 
so as the cj) is increased reflecting the changes in the void 
spaces between particles. To quantify this distribution, 
we obtain the probability distribution function (PDF) of 
A and plot it in Fig. [2ja) as a function of area normalized 
by the median. We find significant changes in the shape 
as a function of 4> with systematically greater deviations 
from F distribution at lower (f). 

To understand these deviations, we note that the par- 
ticle associated with a Voronoi cell appears more asym- 
metrically placed if the particle is located at the edge of 
a large pore. Using this observation, we define a criterion 
based on the distance of the particle from the center of 
the polygonal cell area dy to divide the particles and their 
associated Voronoi areas in to subsets. In Fig. [l] (g-h-i), 
we have indicated a subset of particles with a grey/green 
using a criteria 0.1 < dy/d < 0.5, and a subset using a 
criteria dy/d < 0.1 with red/dark grey. We find that all 
areas with substantial overlap with large void spaces are 
eliminated in this last subset. We then plot PDF of areas 
with dy/d < 0.5 and subset with dy/d < 0.1 along with 
PDF of all the areas in Fig. [2]jb-c-d) for the three rep- 
resentative (j). We find that the PDF for dy/d < 0.1 ap- 
proach a r distribution with ly = 3.6, whereas PDF cor- 
responding to all the cells show strong exponential tails. 
Thus, we are able to link the nature of the distribution of 
free area for all particles to the presence of heterogeneity 
in the system due to the presence of the voids. Addi- 
tional complementary analysis shows that cells become 
more regular shaped with compression and is discussed 
in the supplementary documentation [10]. 

In order to characterize the void spaces or pores di- 
rectly, we investigated the pore size distribution using 
the distance of a random point in a pore to the near- 
est point on the pore-solid interface |16j . As shown in 
the supplementary documentation [10], we find a wide 
distribution of pores size whose average size decreases 
continuously during the compression. This trend reflects 
the evolution of the exponential tail in free area distribu- 
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FIG. 3. (a) The radial correlation function g{r) for ij) = 0.33, 
(j) = 0.50, and — 0.65 (Top to bottom). Experimental re- 
sults are compared with theoretical calculation for hard disks 
given by the Percus-Yevick equation, (b) g{r = d) (amplitude 
of the first peak of this function) versus (j> for experiments with 
attractive particles and for solution of Percus-Yevick equa- 
tion. 



tions ~ exp(— aA), where the parameter a is related to 
the average pore size. 

To examine the short scale structure in the aggregates, 
we evaluate the radial correlation function g[r) which 
measures the particle density-density correlation at dis- 
tance r in Fig.jsja) for various 4>. For hard particles with 
diameter d, g{r) = for r < d, and g{r) — )• 1 for large 
r. Peaks are clearly observed near d, y/id, 2d, \/7d, 3d, 
2y/3d and \/l3d indicating a hexagonal packing. Further 
it may be noted that the locations of the peaks do not 
vary significantly as is increased. These observations 
appear to be consistent with the fact that fioating par- 
ticles are pulled towards each other if they are close by 
until they come in contact. As (j) is increased, the relative 
height of the secondary peaks increase consistent with an 
increase in size of the clusters. 

Now, it is well known that g(r) for non-cohesive hard- 
disks show peaks because of steric interactions as the 
density is increased. Therefore, we also plot in Fig. [sj^a) 
gp.Y. (f) calculated for hard disks arranged randomly us- 
ing the Percus-Yevick equation ^7]. The peak at r — d 
for cohesive particles is significantly higher at low (f) and 
in fact decreases in amplitude with (j) in contrast with 
noncohesive hard disks as also shown in Fig.jsJ^b). It may 
be further noted that the location of the secondary peaks 
does not change with (j), whereas those corresponding to 
the noncohesive disks shift to the left. This is consistent 
with the fact that it is improbable to form a chain of par- 
ticles in contact without cohesion. As cf) increases, g{r) 
appears to grow more similar both in terms of the ampli- 
tude of the first peak and the overall form, indicating the 
increasingly dominant role played by steric interactions 
in organizing the cluster. 

To understand the connectivity of the aggregates, we 
obtain the average number of contacts (z) between par- 
ticles (see Fig. |4][a)). Here, we have used a threshold 
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FIG. 4. (a) Average number of contacts inside capillary ag- 
gregates as a function of (j>- Error bars are estimated from 
statistical dispersion and accuracy of measurement, (b) PDF 
of the local bond order parameter \tp6\, with a peak near 1 
indicating strong hexagonal order. 

distance of separation 1.2d for defining a contact be- 
tween particles to account for cumulative errors in par- 
ticle tracking and to include all particles which fall in 
the first peak in g{r). Particles are isolated in the limit 
^ — 0, and therefore (z) — 0. It may be noted that (z) 
increases rapidly by ^ 0.2 and then does not increase 
significantly over a broad range of (f). Thus, it appears 
that once clusters forms, they only grow larger and their 
local character does not change significantly. The scale 
over which rapid increase in (z) occurs can be estimated 
based on the following arguments. The average Voronoi 
cell is an hexagon of area with side given by the aver- 
age distance between isolated points Ip over V3 [13]. For 
particles randomly distributed inside a circle of radius 

Rq, we get: ttRq'^ = N (Ay) - N^l/ and l/ - 

If Ip < 4c?, because the experiments are conducted slowly 
over 600 s, most of particles are not isolated. Therefore, 
4> ^ 7r/32'\/3 ^ 0.06, (z) becomes of order one, in quali- 
tative agreement with measurements. Finally, disks with 
friction should jam for (z) > 3 according to the Maxwell 
criterion in two dimension |llj . We observe (z) consis- 
tent with this value and well below the value of 6 for a 
hexagonal crystal. 

A complementary and more subtle method to test the 
local order, given the hexagonal order indicated by the 
peaks in g{r), is using the local bond order parameter 

■06 = — ^f=i e*^^'^ . Where, 6ij is the angle between 

the line joining particles i, j and an arbitrary fixed direc- 
tion, and Ni is the number of nearest neighbors for the 
particle i. Nearest neighbors for each particle are identi- 
fied as those within 1.2 d. We then plot the PDF of 
in Fig. [ijjb). A peak is observed at ji/'el = 1 showing the 
prevalence of small hexagonal order. But a broader peak 
is also observed at ^ 0.6 indicating presence of some dis- 
order. Further, the distributions are similar for all the 
low to high (p plotted, also confirming that the packing 
does not evolve significantly with at a local particle 
level. 



In summary, we investigate the aggregation of cohe- 
sive granular matter over a wide range of particle den- 
sity using a novel experimental system and find unan- 
ticipated features in the structures formed. In particu- 
lar, we find that the observed aggregates at low density 
show local close packed structure in contrast with cohe- 
sionless particles, but the overall heterogeneity is signif- 
icantly greater than for cohesionless particles, which we 
characterize with Voronoi cell distributions. Attraction 
is observed to put the system into a local minima of en- 
ergy, which introduces a short range order on a scale of 
few particles. But the structure remains globally disor- 
dered, porous and heterogeneous. By applying uniform 
compression, we show that density is increased by filling 
of pores. Further we find that the local sub-structure of 
the particles is on average conserved during the compres- 
sion. As the density is increased, statistical properties ap- 
proach that for cohesionless frictional granular as steric 
interactions increase in importance. Finally, although 
large scale heterogeneity is reduced at higher densities, 
global order is never reached and the system buckles far 
from an ordered state and at densities below jamming 
threshold for cohesionless particles. 
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